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ABSTRACT 

Type la supernova (SN la), galaxy clustering, and cosmic microwave background anisotropy (CMB) 
data provide complementary constraints on the nature of the dark energy in the universe. We find 
that the three-year Wilkinson Microwave Anisotropy Probe (WMAP) observations give a CMB shift 

parameter of i? = (il,nHQ^ Jn'^''"^ dz' / H{z') = 1.70 ± 0.03. Using this new measured value of the 
CMB shift parameter, together with the baryon acoustic oscillation (BAO) measurement from the Sloan 
Digital Sky Survey (SDSS), and SN la data from the HST/GOODS program and the first year Supernova 
Legacy Survey, we derive model-independent constraints on the dark energy density px{z) and the cosmic 
expansion rate H{z). We also derive constraints on the dark energy equation of state wx{z) — wq + w' z 
(with cutoff at z — 2) and wx{a) = wq + (1 ~ a)wa. 

We find that current data provide slightly tighter constraints on px{z) and H{z) as free functions 
in redshift, and roughly a factor of two improvement in constraining wx{z). A cosmological constant 
remains consistent with data, however, uncertainties remain large for model-independent constraints 
of dark energy. Significant increase in the number of observed SNe la between redshifts of 1 and 2, 
complemented by improved BAO and weak lensing cosmography measurements (as expected from the 
JEDI mission concept for the Joint Dark Energy Mission), will be required to dramatically tighten 
model-independent dark energy constraints. 

Subject headings: Cosmology 



1. INTRODUCTION 

The observed accelerated expansion of the universe 
(Riess et al. 1998; Perlmutter et al. 1999) can be ex- 
plained by an unknown energy component in the universe 
(Freese et al. 1987; Linde 1987; Peebles & Ratra 1988; Wet- 
terich 1988; Frieman et al. 1995; Caldwell, Dave & Stein- 
hardt 1998), or a modification of general relativity (Sahni 
& Habib 1998; Parker & Raval 1999; Dvah, Gabadadze, & 
Porrati 2000; Mersini, Bastero-Gil, & Kanti 2001; Freese 
& Lewis 2002). Padmanabhan (2003) and Peebles & Ra- 
tra (2003) contain reviews of many models. Some recent 
examples of models are presented in Carroll et al. (2004); 
Onemli & Woodard (2004); Cardone et al. (2005); Kolb, 
Matarrese, & Riotto (2005); Martineau & Brandenberger 
(2005); Mclnnes (2005); Cai, Gong, & Wang (2006). For 
convenience, we refer to the cause for the cosmic accelera- 
tion as "dark energy" . 

Data of Type la supernovae (SNe la), cosmic large scale 
structure (LSS), and the cosmic microwave anisotropy 
(CMB) are complementary in precision cosmology (Bah- 
call et al. 1999; Eisenstein, Hu, & Tegmark 1999; Wang, 
Spergel, & Strauss 1999). An important development in 
this complementarity is to use the baryonic acoustic os- 
cillations (BAO) in the galaxy power spectrum as a cos- 
mological standard ruler to probe dark energy (Blake & 
Glazebrook 2003; Seo & Eisenstein 2003). 

In placing robust constraints on dark energy, it is im- 
portant to (1) derive model- independent dark energy con- 
straints (Wang & Garnavich 2001; Tegmark 2002; Daly & 
Djorgovski 2003), and (2) use data derived without assvmv- 
ing dark energy to be a cosmological constant (Wang & 



Tegmark 2004). 

In this paper, we derive the CMB shift parameter from 
the three-year Wilkinson Microwave Anisotropy Probe 
(WMAP) observations (Spergel et al. 2006; Bennett et al. 
2003), and show that its measured value is mostly indepen- 
dent of assumptions made about dark energy. Using this 
new measured value of the CMB shift parameter, together 
with the LSS data from the BAO measurement from the 
Sloan Digital Sky Survey (SDSS), and SN la data from 
the HST/GOODS program and the first year Supernova 
Legacy Survey (SNLS), we derive model- independent con- 
straints on the dark energy density px{z) and the cosmic 
expansion rate H{z). For reference and comparison, we 
also derive constraints on the linear dark energy equation 
of state wx(z) = wq + w'z (with cutoff at z = 2) and 
wx{a) = uiQ -I- (1 — a)wa (Chevallier & Polarski 2001). 

Sec. 2 describes the method and data used in our cal- 
culations. We present results in Sec. 3 and summarize in 
Sec.4. 



2. THE METHOD AND DATA USED 

2.1. The method 

We run a Monte Carlo Markov Chain (MCMC) based 
on the MCMC engine of Lewis & Bridle (2002) to obtain 
O(IO^) samples for each set of results presented in this pa- 
per. The chains are subsequently appropriately thinned. 

We derive constraints on the dark energy density px{z) 
as a free function, with its value at redshifts Zi, px{zi), 
treated as independent parameters estimated from data. 
For z > Zcut, we assume px{z) to be a powerlaw smoothly 
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matched on to px(z) at z 

px{z) = PxiZcut) 



Zcut (Wang & Tegmark 2004) 
1 



z 



(1) 



, 1 + ^^cut , 

The number of observed SNe la is either very few or none 
beyond Zcut- For the Riess et aL (2004) sample, Zcut ~ 1-4. 
For the Astier et aL (2005) sample, Zcut = 1-01. We use cu- 
bic spline interpolation to obtain values of pxiz) at other 
values of z (Wang & Tegmark 2004). 

The H{z) values corresponding to the estimated pxi^i) 
are estimated directly from the MCMC chain to fully in- 
corporate the correlation between the estimated param- 
eters, and compared with the uncorrelated estimates of 
H{z) from SN la data only (Wang & Tegmark 2005). 

For reference and comparison with the work by others, 
we also derive constraints on dark energy models with a 
constant dark energy equation of state w (wxiz) = w), 
and a linear equation of state parametrized by (1) wxiz) = 
Wo -I- w'z at z < 2, and wx{z) = wo + 2w' at z > 2; (2) 
wx{a) = Wo + (1 -a)wa. 

2.2. SN la data 

Calibrated SN la data (Phillips 1993; Riess, Press, & 
Kirshner 1995) give luminosity distances dL{zi) to the red- 
shifts of the SNe la Zi. For a flat universe 

dz' 



dL{z)^cHQ\l + z) 
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E{z 
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il + zy + il-Qm)px{z)/px{ 
with px{z) denoting the dark energy density. 

We use SN la data from the HST/GOODS program 
(Riess et al. 2004) and the first year SNLS (Astier et al. 
2005), together with nearby SN la data. The comparison 
of results from these two data sets provides a consistency 
check. 

We do not combine these two SN la data sets, as they 
have systematic differences in data processing. Fig.l shows 
the difference in estimated distance moduli for 37 SNe la 
included in both the Riess et al. (2004) "gold" sample 
and the Astier et al. (2005) sample; there is clearly signifi- 
cant scatter due to difference in analysis techniques (Wang 
2000b). As a result, the two data sets have noticably dif- 
ferent zero point calibrations. A given zero point calibra- 
tion affects the measuremnt of Hq, but has no impact on 
dark energy constraints (see Eq.[2]). Combining the data 
sets will lead to artificial systematic errors resulting from 
the difference in zero point calibrations, which are diffi- 
cult to quantify; this outweighs the gain in accuracy at 
present since there is a large overlap between the two data 
sets. It is important to use data analyzed using the same 
technique (which corresponds to the same zero point cali- 
bration); although it would be useful to use the same data 
analyzed using one or more other techniques (but only one 
technique should be used at a time) for cross check. 

We use the Riess et al. (2004) "gold" sample fiux- 
averaged with Az = 0.05. This sample includes 9 SNe 



la at z > 1, and appears to have systematic effects that 
would bias the distance estimates somewhat without flux- 
averaging (Wang & Tegmark 2004; Wang 2005). Flux- 
averaging (Wang 2000b; Wang & Mukherjec 2004) removes 
the bias due to weak lensing magnification of SNe la (Kan- 
towski, Vaughan, & Branch 1995; Frieman 1997; Wambs- 
ganss et al. 1997; Holz 1998; Metcaff & Silk 1999; Wang 
1999; Barber et al. 2000; Vale & White 2003)^, or other 
systematic effects that mimics weak lensing qualitatively 
in affecting the observed SN la brightness. 

We have added a conservative estimate of 0.15 mag in 
intrinsic dispersion of SN la peak brightness in quadrature 
to the distance moduli published by Astier et al. (2005), 
instead of using the smaller intrinsic dispersion derived by 
Astier et al. (2005) by requiring a reduced x^ = 1 in their 
model fitting. This is because the intrinsic dispersion in 
SN la peak brightness should be derived from the distri- 
bution of nearby SNe la, or SNe la from the same small 
redshift interval if the distribution in SN la peak bright- 
ness evolves with cosmic time. This distribution is not well 
known at present, but will become better known as more 
SNe la are observed by the nearby SN la factory (Aldering 
et al. 2002) and the SNLS (Astier et al. 2005). By using 
the larger intrinsic dispersion, we allow some margin for 
the uncertainties in SN la peak brightness distribution. 

2.3. LSS data 

For LSS data, we use the measurement of the BAO peak 
in the distribution of SDSS luminous red galaxies (LRG's). 
We do not use the linear growth rate measured by the 2dF 
survey, as there are some inconsistencies in currently pub- 
fished 2dF resuhs (Verde et al. 2002; Hawkins et al. 2003).^ 

The SDSS BAO measurement (Eiscnstein et al. 2005) 
gives A = 0.469 {ns/ 0.98)-"-^^ ± 0.017 (independent of a 
dark energy model) at zbao — 0.35, where A is defined as 
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r^{zBAo) 



czbao 
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CZBAO 
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ii{zBAo) _ 

where r(z) is the comoving distance, and -ff (z) is the Hub- 
ble parameter. Note that ii{z) = HoE{z). We take the 
scalar spectral index ns = 0.95 as measured by WMAP3 
(Spergel et al. 2006). Note that this constraint from Eiscn- 
stein et al. (2005) is not just a simple measurement of the 
BAO feature; it also relies on the constraints on fimh^ 
from measuring the power spectrum turnover scale (re- 
lated to matter-radiation equality) . The latter makes their 
BAO constraint less robust than it would be otherwise. A 
new analysis of the SDSS data to derive truly robust BAO 
constraints would be very useful for placing dark energy 
constraints (Dick, Knox, & Chu 2006). 

Also note that the Eiscnstein et al. (2005) constraint on 
A depends on the scalar spectral index ns- Since the error 
on ns from WMAP data does not increase the effective er- 
ror on A, and the correlation between ns and the CMB 
shift parameter R is weak, we have ignored the very weak 
correlation between A and R in our likelihood analysis. 

sec Dodclson & Vallinotto (2005); Cooray, Holz, & Huterer 



^ Weak lensing magnification of SNe la can also be used as a eosmological probe, 
(2006); Munshi & Valageas (2006). 

■* Combining the results from Verde et al. (2002) and Hawkins et al. (2003), we would obtain the linear growth rate f{z2df) = 0.51 ± 0.11. 
However, Hawkins et al. (2003) points out that the Verde et al. (2002) results strongly depend on the assumed pairwise peculiar velocity 
dispersion of 385 km/s, while Hawkins et al. (2003) finds the pairwise peculiar velocity dispersion to be 500 km/s. A more self-consistent linear 
growth rate from 2dF data has to await a new bispeetrum analysis of the 2dF data. 
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We have derived R from WMAP data marginalized over 
all relevant parameters. 

2.4. CMB data 

The CMB shift parameter R is perhaps the least model- 
dependent parameter that can be extracted from CMB 
data, since it is independent of Hq. The shift parameter 
R is given by (Bond, Efstathiou, & Tegmark 1997) 



i? = fi™'/' 



dz'/E{z'), 



(5) 



where zcmb is the redshift of recombination; thus R = 
(jlmHQ) r(zcMB)/c (for a flat universe), which is well 
determined since both flmh^ and t^zcmb) are accurately 
determined by CMB data. Similar reasoning applies to an 
open and a closed universe as well. 

The ratio of the sound horizon at recombination to 
r{zcMB), 6s, should be more or less equivalent to R. We 
call R the "CMB shift parameter" following previously 
published literature over the last 9 years (see for exam- 
ple. Bond, Efstathiou, & Tegmark (1997); Odman et al. 
(2003)). One advantage of R is that it only involves a sim- 
ple integral over 1/E{z), while the sound horizon at last 
scattering is more complicated to calculate accurately and 
depends on more parameters. 

We compute R using the MCMC chains from the analy- 
sis of the three year WMAP data provided by the WMAP 
team (Spergel et al. 2006). The resultant probability dis- 
tribution of R for three different classes of models are 
shown in Fig. 2. Clearly, the measured R from WMAP 
3 year data has only a very weak model dependence^ . 

We use a Gaussian distribution in R with (R) — 1.70 
and an = 0.03 (thin solid line in Fig. 2) in deriving our re- 
sults presented in Sec. 3. This fits the WMAP 3 year data 
well, and allows some margin in error for the very weak 
model dependence of R. 

3. RESULTS 

In deriving all the results presented in this section, we 
have assumed a flat universe, and marginalized over Q„i 
and h. 

Fig. 3 shows px{z) measured using SN la data (Riess et 
al. 2004; Astier et al. 2005), combined with the WMAP 
3 year data, and the SDSS BAO data. Note that beyond 
•2^cut=l-4 (upper panel) and 1.01 (lower panel), px{z) is 
parametrized by a power law (1 -I- z)", the index of which 
is marginalized over. 

Fig. 4 shows the uncorrelated H{z) estimated using only 
SN la data (Wang & Tegmark 2005), and the H{z) corre- 
sponding to the px{zi) shown in Fig. 3. 

Tables 1 and 2 give the 68% confidence intervals of 
Px{zi), a, and H{zi) for the WMAP 3 year and the 
SDSS BAO data, combined with SN la data from Riess 
et al. (2004) and Astier et al. (2005) respectively. 

Tables 3 and 4 give the covariance matrices for (px{zi), 
a), and {H{z,), a) for the WMAP3 and SDSS BAO data 
combined with SN la data from Riess et al. (2004) and 
Astier et al. (2005) respectively. 

Table 5 gives the constraints on a constant w 
{wx{z) =const.), {wo,w') {wx{z) = wq + w'z at z < 2, 



and wx{z) ^ wq + 2w' at z > 2), and (wq, Wa) (wxia) = 
Wo + (1 - a)wa). 

Fig. 5 shows the 68% and 95% joint confidence contours 
for {wo,w') and (wo^Wa)- 



4. DISCUSSION AND SUMMARY 

In order to place robust constraints on dark energy 
in a simple and transparent manner, we have derived 
the CMB shift parameter from the WMAP 3 year data 
{R = 1.70 ± 0.03, see Fig. 2). We constrain dark energy 
using this new measurement of the CMB shift parameter, 
together with LSS data (the BAO measurement from the 
SDSS LRG's), and SN la data (from the HST/GOODS 
program and the first year SNLS). 

We have derived model-independent constraints on the 
dark energy density px{z) and the cosmic expansion rate 
H(z) (Figs. 3-4 and Tables 1-4). There are two reasons 
that one should use pxiz) and H{z) instead of wx{z) 
to probe dark energy. First, px{z) and H{z) are more 
directly related to observables than wxiz) (see Eqs.[2]- 
[5], and note that d\npx/d\na = —3(1 -I- wx)) (Wang 
& Garnavich 2001; Tegmark 2002). This means that 
px{z) is more tightly constrained by data than wx{z) 
(Wang & Freese 2006; Daly & Djorgovski 2005; Huterer 
& Cooray 2005; Ishak 2005; Dick, Knox, & Chu 2006). 
Secondly, px{z) and H{z) are more general phenomeno- 
logical representations of dark energy than wx{z) (Wang 
& Tegmark 2004). One must integrate the equation, 
d\npx/dhia =■ —3(1 -I- wx\ to obtain pxiz) before a 
comparison with data can be made, hence even arbitrary 
function wxiz) has the hidden assumption that px{z) is 
non- negative (since /9x(0) > 0). Since we don't know what 
dark energy is - it may not even be energy at all but a 
modification of general relativity - px{z) may well have 
been negative at some past epoch or become negative in 
the future. Measuring px{z) and H{z) (instead oi wx{z)) 
from data allows us to constrain this possibility. 

For comparison with the work by others, we also derive 
constraints on the dark energy equation of state wx (z) = 
wq+w'z (with cutoff at z = 2) and wx{a) = wo + {i — a)wa 
(Fig.5 and Table 5). 

Because of the difference in SN la analysis techniques 
used in deriving the Riess et al. (2004) "gold" sample and 
Astier et al. (2005) sample (see Fig.l), we have presented 
the results for these two data sets separately. This is nec- 
essary in order to avoid introducing systematic errors due 
to the different zero point calibrations of the two data sets 
(which are difficult to quantify) (Wang 2000b). We find 
that the Riess et al. (2004) and Astier et al. (2005) data 
sets give similar and consistent constraints on dark energy 
(Figs. 3-5), with the Riess et al. (2004) sample being able 
to constrain px{z) and H{z) at z > 1 since it contains 9 
SNc la at z > 1. This is reassuring, and seems to indicate 
that the Riess et al. (2004) sample (after flux-averaging) 
is consistent with the Astier et al. (2005) sample. It will 
be useful to derive the distances to all the SNe la using the 
same technique to reduce systematic uncertainties (Wang 
2000b) and maximize the number of SNe la that can be 
used in the same dark energy analysis. 



^ For a Bayesian analysis of the number of parameters required by current cosmological data, in other words of what comprises an adequate 
model, see Mukherjee, Parkinson, & Liddle (2006) 
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Table 1 

The mean and the 68% confidence intervals of px(zi), a, and H{zi) for the CMB and LSS data combined with SN 

Ia data from Riess et al. (2004). 



Parameter 


Riess04+WMAP3+SDSS 


px (0.467) 

px (0.933) 

Px (1.400) 

a 


1.159 (0.953, 1.361) 

1.357 (0.602, 2.095) 

2.751 (0.408, 5.053) 

-0.037 (-1.927, 1.818) 


J?(0.467) 
i7(G.933) 
iJ(1.400) 


1.321 (1.257, 1.384) 
1.758 (1.617, 1.896) 
2.419 (2.045, 2.797) 



Table 2 

The mean and the 68% confidence intervals of px{zi), a, and H(zi) for the CMB and LSS data combined with SN 

Ia data from Astier et al. (2005). 



Parameter 


Asticr05+WMAP3+SDSS 


px(0.505) 

Px(l.OlO) 

a 


1.013 (0.893, 1.131) 
1.579 (0.833, 2.327) 
0.030 (-1.786, 1.770) 


ff(0.505) 
ff(l.OlO) 


1.290 (1.252, 1.327) 
1.828 (1.678, 1.978) 



Table 3 

The covariance matrices for px{zi) and H{zi) from the WMAP 3 year and the SDSS BAO data combined with 

SN Ia data from Riess et al. (2004). 





Px (0.467) 


px (0.933) 


Px (1.400) 


a 


pjf (0.467) 

px (0.933) 

Px (1.400) 

a 


0.426E-01 
0.226E-01 
-0.137E+00 
0.153E-01 


0.226E-01 
0.585E+00 
-0.261E+00 
-0.322E+00 


-0.137E+00 
-0.261E+00 
0.717E+01 
-O.llOE+01 


0.153E-01 
-0.322E+00 
-O.llOE+01 
0.276E+01 




i?(0.467) 


iJ(0.933) 


H{1A00) 


a 


i?(0.467) 

iJ(0.933) 

H{1A00) 

a 


0.407E-02 
0.185E-02 
0.121E-02 
0.851E-02 


0.185E-02 
0.200E-01 
-0.130E-02 
-0.534E-01 


0.121E-02 
-0.130E-02 
0.154E+00 
-0.125E+00 


0.851E-02 
-0.534E-01 
-0.125E+00 
0.276E+01 



We find that compared to previous results in Wang & 
Tegmark (2004), the current data provide a slightly tighter 
constraint on px{z) and H{z) as free functions of cosmic 
time, and roughly a factor of two improvement in con- 
straining wx{z)- Note that the Astier et al. (2005) data 
set (together with WMAP 3 year and SDSS BAO data) 
provides a tighter constraint on wq = wxiz = 0) because 
it contains more SNe la at lower redshifts. Because of the 
strong correlation between w' (or Wa) and wq, this has led 
to tighter constraints on w' and Wa as well. 

A cosmological constant remains consistent with data, 
although more exotic possibilities are still allowed, con- 
sistent with previous results (see for example, Wang 
& Tegmark (2004, 2005); Alam & Sahni (2005); Daly 



& Djorgovski (2005); Jassal, Bagla, & Padmanabhan 
(2005a,b); Dick, Knox, & Chu (2006); Ichikawa & Taka- 
hashi (2006); Jassal, Bagla, & Padmanabhan (2006); 
Nesseris & Perivolaropoulos (2006); Schimd et al. (2006); 
Wilson, Chen, & Ratra (2006)). In particular, Zhao et 
al. (2006) uses WMAP 3 year, SN la data from Riess et 
al. (2004), together with SDSS 3D power spectra and 
Lyman-a forest information data, and found constraints 
on wx{z) = wo + wiz/{l + z) that are qualitatively consis- 
tent with our results, with significant differences that can 
be explained by the differences in the combination of data 
used, and perhaps some data analysis details. 

We find that uncertainties remain large for model- 
independent constraints of dark energy (see Figs. 3-4). Sig- 
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Table 4 

The covariance matrices for px{zi) and H(zi) from the WMAP 3 year and the SDSS BAO data combined with 

SN Ia data from Astier et al. (2005). 





px(0.505) /9x(l-010) a 


px(0.505) 

Pjf(l.OlO) 

a 


0.148E-01 0.472E-01 -0.750E-01 
0.472E-01 0.559E+00 -0.460E+00 
-0.750E-01 -0.460E+00 0.254E+01 




H{0.505) iJ(l.OlO) a 


7?(0.505) 

7?(1.010) 

a 


0.146E-02 0.308E-02 -0.118E-01 
0.308E-02 0.223E-01 -0.729E-01 
-0.118E-01 -0.729E-01 0.254E+01 



Table 5 
The mean and 68% and 95% confidence level constraints on a constant w and {wi),w'), and the covariance 

between Wo and w' , and between KJo and Wa- 





Ricss04+WMAP3+SDSS 


Asticr05+WMAP3+SDSS 


w 


n SSCi +0.109 +0.206 
"•°°^ -0.111 -0.227 


r, qqq +0.082 +0.159 
u.yyy -o.083 -0.I68 


Wo 


(. 704 +0.243 +0.584 
-U./yi -0.244 -0.431 


n QSQ +0.100 +0.413 

-u.yoy -0.162 -0.291 


w' 


_n 446 +0.711 +0.948 
U.11U _o.710 -2.237 


n 1 77 +0.571 +0.742 
"•-^' ' -0.574 -1.984 


(j'^{wo,w') 


-0.192 


-0.112 


Wo 


n 01 +0.293 +0.704 
U.OiO _o.296 -0.508 


1 ni7 +0.199+0.503 
i.UK -0.200-0.350 


Wa 


^Sin +i-^05 +1.792 
-U.OiU _i.259 -3.620 


n nao +1.045 +1.429 

-U.UOy -1.052 -3.173 


a^{wo,Wa) 


-0.410 


-0.245 



nificant increase in the number of observed SNe la between 
redshifts of 1 and 2, complemented by BAO and weak lens- 
ing cosmography measurements (such as expected from 
the JEDl mission concept of the Joint Dark Energy Mis- 
sion), should dramatically tighten dark energy constraints 
and shed light on the nature of dark energy (Wang 2000a; 
Wang et al. 2004). ^ 
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Fig. 1. — The difference in estimated distance moduli for 37 SNe la included by both the Riess et al. (2004) "gold" sample and Astier et 
al. (2005); there is clearly significant scatter due to difference in analysis techniques (Wang 2000b). 
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